Topological entropy of realistic quantum Hall wave functions 
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The entanglement entropy of the incompressible states of a realistic quantum Hall system are 
studied by direct diagonalization. The subdominant term to the area law, the topological entangle- 
ment entropy, which is believed to carry information about topologic order in the ground state, was 
extracted for filling factors f = 1/3, v — 1/5 and v = 5/2. The results for v — 1/3 and u — 1/5 are 
consistent with the topological entanglement entropy for the Laughlin wave function. The v = 5/2 
state exhibits a topological entanglement entropy consistent with the Moore-Read wave function. 
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I. INTRODUCTION 

This paper is a numerical study, using direct diagonal- 
ization, of the entanglement entropy of incompressible 
states of quantum Hall systems. In particular, the en- 
tanglement entropy is calculated for filling factors 1/3 
and 1/5 in the n = Landau level and the 5/2 state in 
the n = 1 Landau level. The primary motivation for this 
work is to better understand the nature of the 5/2 state 
observed in experimenlji. Originally, the 5/2 state was 
not believed to be spin polarized^, however, with theo- 
retical inpuli^i^i^ and further experimental investigation^ 
an incompressible spin polarized state was revealed. An 
elegant theoretical possibility for this state is the wave 
function suggested by Moore and Read^; the novel fea- 
ture of the Moore-Read state being the presence of non 
abelian fractional statistics. As well as being of inter- 
est in its own right, states with non abelian statistics 
give rise to possible robust implementation of quantum 
computation^'*. However, the Moore- Read state is the 
ground state of a not very realistic Hamiltonian with a 
three body interaction term^. 

The question remains, does the Moore-Read state con- 
tain the physics of the v — b/2 quantum Hall system; 
that is, is it in some sense close to the ground state of 
a model with a realistic Hamiltonian, i.e. electrons in 
a magnetic field interacting via long range Coulomb in- 
teraction? A valuable way to address this issue is to 
compare the ground state wave function obtained from 
direct diagonalization of a realistic Hamiltonian to the 
Moore-Read wave function. The results of such studies 
are unfortunately somewhat ambiguous^iiiS. In addition, 
even if there is a large overlap between the numerical 
wave function and the Moore-Read state, how does one 
know whether this overlap truly indicates that the long 
distance, low energy behavior is the same, in particu- 
lar whether the numerical wave function has non abelian 
statistics. Attempts to directly detect non abelian statis- 
tics in numerical systems with strictly Coulomb interac- 
tions have thus far also proven to be elusiv o^^'^^ . Note 
in reference [12], non abelian statistics were clearly ob- 
served, however, the Hamiltonian contained a mixture of 
Coulomb and 3-body terms. 



Recently, a novel approach, using concepts from quan- 
tum information theory, has been proposed to character- 
ize incompressible quantum states. Kitaev and PreskilU^ 
and Levin and WeniA have shown the sub leading contri- 
bution to the entanglement entropy of a subsystem, the 
topological entanglement entropy is universal and reflects 
the statistics of quasiparticles of the incompressible Hall 
state in question. To be more explicit, consider a two 
dimensional quantum many body system and spatially 
divide the system into two parts, the part of interest be- 
ing the subsystem and the rest which is referred to as 
the environment. The result of ref. [13,14] is that the 
entanglement entropy scales as 



S~aL--i + 0{\ 



(1) 



for topologically ordered states. Here L is the linear size 
of the boundary of the subsystem and 7 is the topolog- 
ical entanglement entropy. S, the entanglement entropy, 
more precisely defined in section II, intuitively is a mea- 
sure of the quantum entanglement of the subsystem and 
the environment. The sub leading term, the topological 
entanglement entropy was showni^ii^ to be equal to the 
log of the total quantum dimension of the state in ques- 
tion. Intuitivelyfi^ 7 reflects the number of distinct quasi 
particle types and how the number of linear independent 
states, for each type, grows with the number of quasi 
particles in the states. 

It is important to note that the above expression holds 
asymptotically for large L, that is, a large subsystem and 
a very large environment. By considering several differ- 
ent subsystems, it was shown that the leading contribu- 
tion, which scales as the linear size of the subsystem (and 
is non universal) could be cancelled out and the sub lead- 
ing term could be extracted. Numerically, however, it is 
not easy to implement this methodic. 

In references [17] and [18], a more practical method for 
quantum Hall systems, based on looking at several sys- 
tem sizes, was used to successfully compute the topolog- 
ical entanglement entropy of the Laughlin wave function 
for 1/ — 1/3 and the Moore-Read state for = 5/2. In 
this paper we apply the method of reference [17] to the 
exact ground state wave functions obtained from realis- 
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tic Hamiltonians by direct diagonalization. If the value 
of the topological entanglement entropy extracted from 
a direct diagonalization calculation agrees with the value 
calculated for the Moore-Read state this provides evi- 
dence that the Moore-Read state correctly describes the 
physics of filling factor v ~ b/2. It is important to re- 
alize, without other physical constraints, equality of the 
topological entanglement entropy is a necessary but not 
sufficient condition for states to be topologically equiv- 
alent. Simply stated, if two states don't have the same 
value of 7, they are not equivalent. However, equality of 
7 does not necessarily imply two states are equivalent. 
Therefore to completely characterize an incompressible 
state additional information is needed. 

The paper is organized in the following way: In the 
next section, the numerical method is described, high- 
lighting the differences from reference [17], and in the fol- 
lowing section the numerical results are presented. The 
final section is a summary and gives our conclusions. 



II. NUMERICAL METHOD 

To do our direct diagonalization (DD) calculations, we 
work with finite square clusters with periodic boundary 
conditions, the flat torus geometry. This is in contrast 
to reference [17] which works in the spherical geometry. 
One reason we favor the torus geometry is that most 
published density matrix renormalization group (dmrg) 
calculations of quantum Hall systems, which can handle 
larger system sizes, are performed in this geometryi2,(see 
however, Feguin et ali^ for dmrg in the spherical geome- 
try) . Although the present work is strictly DD, we hope 
to lay the ground work for a future dmrg study. 

To take the magnetic field into account, the Landau 
gauge is chosen, where the momentum in the y direction 
is a conserved quantity and the single particle orbitals, 
strips of width the magnetic length, are oriented parallel 
to the y-axis. The single particle orbitals are labeled 
by the x guiding center coordinate or equivalently the 
momentum in the y direction. Although momentum in 
the X direction is also conserved, for simplicity, and since 
this symmetry has not been implemented in dmrg, we do 
not make use of this quantum number—. The Lanczos 
algorithm is used to calculate the ground state in each 
sector of total y momentum and then the lowest energy 
state is selected. For the v = \/Z and v — 1/b fillings 
the state space is restricted to the n = Landau level 
while for the j/ = 5/2 filling the state space is restricted 
to the n = 1 Landau level. We again emphasize that the 
electrons interact via the long range Coulomb interaction 
periodically continued in the usual way^^. Due to cpu 
limitations the calculations were limited to system sizes 
smaller than or equal to 12 electrons in 36 orbitals {ly — 
1/3), 8 electrons in 40 orbitals {v = 1/5), and 16 electrons 
in 32 orbitals, — 5/2). ( the diagonalizations involve 
state spaces of sizes at most 35 x 10^. It is possible, with 
some difficulty to extend the calculations by 1 electron 



for p = 1/5 ). 

We now turn to the method used to compute the 
topological entanglement entropy. After computing the 
ground state (or states, the issue of ground state de- 
generacy will be addressed later) the entanglement en- 
tropy is then calculated. As in the spherical geometry, in 
the torus geometry, there is a natural (numerically easy) 
choice for the subsystem to calculate the entanglement 
entropy. The subsystem chosen consists of / adjacent (in 
x) orbitals, for example, for / = 2 one can take orbitals 
1 and 2 to get the 2 orbital entanglement entropy. In a 
system with N total orbitals (i.e. for 11 electrons in 33 
orbitals, N — 33) the many electron wave function has 
the form of a collection of coefficients ^iii2i3...i„ where 
ii, Z2, is ■ ■ • *7V take the values or 1 and 4" is the ampli- 
tude for the state with occupancies ii,i2,iz ■ ■ - in ■ One 
then computes the I orbital density matrix (an object 
very familiar from the density matrix renormalization 
group^ ). Exphcitly for 1 — 2, 

The density matrix is then diagonalized, yielding the 
eigenvalue Xj from which the entanglement entropy 

-^A.lnA, (3) 
i 

is obtained. 

Since the torus geometry is used in calculating the, 
say, 2 particle entropy, it does not matter if one takes, 
the orbitals 1,2 or 30, 31 etc. However, there is a sub- 
tlety in that we are not explicitly taking conservation of 
x-momentum into account and due to ground state de- 
generacy the numerical wave function may not be trans- 
lationally invariant. To handle this problem, we have 
used the wave functions that are more translationally in- 
variant in the following sense: (take for concreteness 7 
electrons in 21 states, = 1/3 ). Calculate the "left" 
and "right" / body entanglement entropies, taking the 
"left" subsystem to be orbitals 1, 2, 3, . . . , Z and the right 
subsystem to be 22 — Z, . . . , 20, 21. We choose the ground 
state where the entanglement entropies for the " left" and 
"right" subsystems are equal. For example, for 7 elec- 
trons in 21 orbitals, the states with ky = 7, 14, 21 (in ap- 
propriate units) are degenerate but only ky = 14 satisfies 
the above criteria. For = 1/3 and v = 1/5, this criteria 
uniquely picks the ground state. For v = 5/2, there are 
certain filling factors (i.e. 14/28) where 2 ground states 
satisfy the equality of the "left" and "right" entangle- 
ment entropies. In these cases, we pick the ground state 
with the lowest momentum. 

Following reference [17], from the entanglement en- 
tropies for different numbers of orbitals and different sys- 
tem sizes we have attempted to calculate the topological 
entanglement entropy. The idea is based on the asymp- 
totic formula eq.(l). The first term in eq. (1) is re- 
ferred to as the area law, entropy being proportional to 
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FIG. 1: Entanglement entropy of v — 1/3 state, (a) Entan- 
glement entropy versus a/I (where / is number of orbitals com- 
prising the sub system) for different total number of orbitals, 
N = 21- 36. Finite size scaling (ESS), N ^ oo, results for 
S versus v7 (•) is shown with linear least squares fit (LSF), 
yielding an intercept of —7 = —1.13 ± 0.38. (b) Finite size 
scaling of entanglement entropy. N 00 results, shown in 
figure la, were obtained by linear LSF to N = 24, 27, 30, 33, 36 
data. Estimated uncertainty {as) in A?" — > 00 values of S was 
as < 0.1, smaller than the size of the plotting symbol (•). 



least the largest system we can compute with DD) and a 
sufficiently large but not too large subsystem. As in ref. 
[17], we identify the number of Landau orbitals with the 
area enclosed. 
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FIG. 2: Entanglement entropy of = 1/5 state. Same as 
figure 1 except; (a) N ^ 00 results for S versus \/l (•) yield 
an intercept of —7 — —1.62 ± 0.16 (b) N ^ 00 results were 
obtained by linear LSF to N = 25, 30, 35, 40. Estimated un- 
certainty in iV — » 00 values of S was as < 0.15, approximately 
the size of the plotting symbol (•). 



the bounding area, which is one-dimensional in this case. 

The system and subsystem miist be sufficiently large to 
realize the bulk behavior of the strongly correlated state, 
and the subsystem must be much smaller than the sys- 
tem to realize the area law behavior of the entanglement 
entropy. Consider first of all a very large system (or at 



By plotting S vs the square root of the number of or- 
bitals, one hopes to get a straight line (reflecting the 
area law eq. (1)); then the y-intercept should give mi- 
nus one times the topological entanglement entropy. It 
is however, not easy to treat a very large system with di- 
rect diagonalization. A possibility to overcome this dif- 
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ficulty, is to use information from a number of system 
sizes. Again following [17] one can plot the I orbital en- 
tanglement entropy vs {N = number of orbitals) and 
try to extrapolate to an infinite system. Our experience 
indicates that a linear extrapolation in 1 /N works better 
than adding nonlinear terms, so we use a simple linear 
extrapolation only. The extrapolated Z-orbital entangle- 
ment entropies are then plotted vs the square root of the 
number of orbitals in the subsystem (v^). If a straight 
line results, then the topological entanglement entropy 
can be extracted as the y-intercept. 

A possible alternative approach to the above method is 
to numerically implement the method of ref. [13] or ref. 
[14]. Such an approach was taken in ref. [16] and applied 
to the quantum dimer model. Since the dimer model is 
defined on a lattice, it is easier to vary the choice of the 
subsystem. However, to obtain accurate results, large 
system sizes had to be used (this is practical since no 
diagonalization is needed at the Rokshar-Kivelson point 
of the phase diagram) or alternatively special properties 
of the dimer model had to be utilized. 



III. NUMERICAL RESULTS 

Let us now examine the results of our calculation. First 
consider i' — 1/3 in figures la, lb. Figure la is a graph 
of 5, the entanglement entropy, vs the square root of the 
number of orbitals in the subsystem, for various system 
sizes. We have also shown in this figure an extrapola- 
tion to large N given by the solid circles. (See below 
for discussion). For the largest system N = 36, one sees 
reasonably linear behavior for up to 6 orbitals in the sub- 
system and a y-intercept of w — 1 and a topological en- 
tanglement entropy of « 1. In the spherical geometry, 
for ly = 1/3, it is known that the topological entangle- 
ment entropy takes a value of 7 ~ 0.55. However, for the 
torus, the subsystem has two boundaries rather than one 
and even the leading term in the entanglement entropy 
only depends on orbitals rather close (a/I) to the respec- 
tive boundaries. Hence one expects the contribution from 
each boundary to merely add"'^'* and give twice the value 
of topological entanglement entropy for the sphere. Our 
numerical value 7 « 1, given the uncertainties, see figure 
la, is consistent with this expectation. 

To try to get a more precise estimate in figure lb we 
plot S, the Z-orbital entanglement entropy, vs 1/A^ for 
various / values. By doing a linear least squares fit in 1 /N 
for N=24,27,30,33,36 and taking the value of the line at 
1 /N = the Z-orbital entanglement entropy was extrap- 
olated to large system sizes. These values are shown by 
the solid circles in figure la. The solid circles are then 
fit to a straight line and the intercept gives a topological 
entanglement entropy 7 = 1.13 ± 0.38. This is consis- 
tent with the value expected from the Laughlin state for 
1/ — 1/3. The large uncertainty is the result of the poor 
linearity of a curve passing through the solid circles. 

In figure 2a, b we plot analogous graphs for v = 1/5. 



Looking at figure 2a, the entanglement entropy, vs the 
square root of the number of orbitals in the subsystem, 
for the largest system size {N = 40, but only 8 electrons) 
we again see linear behavior for up to 6-7 orbitals. A 
fit to the linear part gives a y-intercept of about -1 or 
a slightly smaller value -1.2 if we exclude 1 = 1. This is 
compared to an expected value of —2 In \/5 — —1-61. The 
solid circles in figure 2a are values of the entanglement 
entropy extrapolated from figure 2b. A linear fit to the 
extrapolated values yields an intercept of —1.62 ± 0.16. 
This is in excellent agreement with the value from the 
I' = 1/5 Laughlin state. The much smaller, though still 
substantial, uncertainty reflects the better linearity of a 
curve passing through the solid circles in comparison to 
the case i/ = 1/3. 

Finally, let us investigate the z/ = 5/2 systemiS. In fig- 
ure 3a we again plot , the entanglement entropy, vs the 
square root of the number of orbitals in the subsystem. 
There is reasonably linear behavior up to about 6-7 or- 
bitals yielding a y-intercept of roughly -1.5. On the other 
hand, for the Moore-Read state we expect a topological 
entanglement entropy of about 21nV8 ~ 2.08. To try 
to get a more precise estimate in figure 3b we plot the l- 
orbital entanglement entropy vs 1 /N for various values of 
I. Again a simple linear extrapolation of the five largest 
systems sizes is adequate. The solid circles in figure 3a 
shows the extrapolated entanglement entropies. Using a 
linear fit, the extrapolated values then give a y-intercept 
of —2.01 ± 0.19. Thus, the numerical state has a topo- 
logical entropy close to that of the Moore-Read state. 
A possible alternative state, having the same topological 
entanglement entropy, is the Halperin 3-3-1 state (We 
thank one of the referees for pointing this out to us). 
However, since we are only working with spin polarized 
electrons in a single layer, this state is excluded on gen- 
eral grounds of symmetry^^ . 



IV. CONCLUSION 

In this paper, direct diagonalization, by necessity on 
small system sizes, has been used to calculate the entan- 
glement entropy for the v = 1/3, v = 1/5 and = 5/2 
quantum Hall states. We emphasize that a realistic 
Hamiltonian, long range Coulomb interaction, has been 
used. At all filling fractions considered, the area law 
has been verified by examining the largest system size. 
To accurately extrapolate the topological entanglement 
entropy it was necessary to extrapolate the I orbital en- 
tanglement entropy to large system sizes. For = 1/3 
and v = 1/5, this extrapolation gave results consistent 
with values of the topological entanglement entropy for 
the Laughlin state. It should be noted that our results 
are consistent with an entropy corresponding to a topo- 
logical ground state with two boundaries (following from 
our computation in the Landau gauge) rather than a sin- 
gle boundary in the spherical geometry. For v = 5/2 
the value of the topological term obtained was consistent 
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with the topological term for the Moore-Read state. We 




view this as a confirmation, that the incompressible state 
a,t u = 5/2, is in fact, the Moore- Read state. 
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FIG. 3: Entanglement entropy oi p = 5/2 state. Same as 
figure 1 except: {a) N oo results for S versus a/I (•) yield 
an intercept of —7 = —2.01 ± 0.19 (b) N ^ 00 results were 
obtained by linear LSF to N = 24, 26, 28, 30, 32. Estimated 
uncertainty in iV — » 00 values of S was as < 0.18, approxi- 
mately the size of the plotting symbol (•). 
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